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Abstract 

We formulate the thermodynamic Bethe Ansatz (TBA) equations for the closed (pe- 

(2) 

riodic boundary conditions) A2 quantum spin chain in an external magnetic field, in the 
(noncritical) regime where the anisotropy parameter rj is real. In the limit 77 — > 0, we 
recover the TBA equations of the antiferromagnetic sM(3)-invariant chain in the funda- 
mental representation. We solve these equations for low temperature and small field, and 
calculate the specific heat and magnetic susceptibility. 
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1. Introduction and Summary 



Given a quantum integrable lattice model in one space dimension, one can find the 
eigenvalues of the Hamiltonian in terms of solutions of the model's Bethe Ansatz (BA) 
equations. (See, e.g., Ref. 1.) Unfortunately, having found the eigenvalues, one is still 
quite far from determining the model's physical properties. The main reason for this 
is that the BA equations are in general very difficult to solve, in particular for N (the 
number of lattice sites) finite. Considerable simplification occurs in the thermodynamic 
{N — > oo) limit. Provided one can formulate a suitable "string hypothesis" for the solutions 
of the BA equations, the problem is then to determine the densities Pn{^) and Pn{^) of 
quasi-particles and quasi-holes, respectively. In principle, this can be accomplished once 
one solves the so-called thermodynamic Bethe Ansatz (TBA) equations for the quantities 
e„ = T In (pn/ Pn)- Since the TBA equations are an infinite set of coupled nonlinear integral 
equations, in practice one solves them perturbatively (e.g., near T — 0). 

This program^ has been successfully applied to a number of integrable lattice models. 
Foremost among these are integrable quantum spin chains - e.g., the spin 1/2 Heisenberg 
chain and its many generalizations. The large body of work^~^^ on the thermodynamics of 
quantum spin chains has had significant consequences for both quantum field theory and 
condensed matter physics. (For a recent introduction, see Ref. 14.) 

In this paper, we focus on the closed (i.e., periodic boundary conditions) chain, 
with Hamiltonian 



where V is the permutation matrix, and R{u) is the i?-matrix associated with the twisted 
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N-l 





(1.1) 



Here 



R{u) = VR{u) , 



(1.2) 



affine algebra A2 in the fundamental representation, which depends on the so-called 
anisotropy parameter rj. The Hilbert space is ®^V, where V is three dimensional. 

This model was first explicitly constructed by Izergin and Korepin^^. The spectrum 
of the transfer matrix and the BA equations were first determined using the analytical BA 
method^^, and later using the algebraic BA method^^. The corresponding vertex model is 
equivalent to an 0{n) model on a square^® or hexagonal^^ lattice. Such 0{n) models are 
relevant^^ to the study of polymers. 

We distinguish two regimes: rj purely imaginary and 77 purely real. For zero external 
magnetic field, these correspond to critical and noncritical regimes, respectively. By abuse 
of language, we shall refer to the two regimes as "critical" and "noncritical" even for 
nonzero field. 

(2) 

For the A2 chain in the critical regime, a general string hypothesis has not yet been 
formulated. A new 2-string solution was found in Ref. 21, and several new candidate 
4-string solutions were found in Ref. 22. Presumably, there are new longer strings as 
well. In the absence of a suitable string hypothesis, the TBA equations of the critical A2 
chain cannot be formulated. An alternative approach of investigating this model, based 
on finite-size corrections, has been recently pursued by two groups^"^'^^. However, there is 
some disagreement between their results. 

The difficulty in formulating a string hypothesis may be related to the fact that, 

(2) 

in the critical regime, the Hamiltonian of the A2 chain is not Hermitian. One might 
try to restrict the space of states to a subspace in which the Hamiltonian is Hermitian. 
However, since this model does not have a quantum-algebra symmetry, it is not clear how 
to implement this restriction. 

The situation for the open A^"* chain, whose Hamiltonian is the same as (1.1) except 
without the final term TCn,!, may be better. As shown in Ref. 25, this model has the 
quantum algebra symmetry Uq[su{2)], and is integrable. Indeed, the spectrum of the 
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transfer matrix and the BA equations have been found, using a generahzation of the 
analytical BA method, in Ref. 26. In the critical regime, one may be able to exploit the 
model's Uq[su{2)] symmetry to make suitable projections on the space of states, in analogy 
with the A^^^ case^'''. 

As a warm-up exercise before addressing these problems, we consider in this paper 
the closed A2 chain in the noncritical* regime. Here the Hamiltonian is Hermitian, and 
there is no difficulty in formulating^^ a string hypothesis. We probe this system at finite 
temperature T in an external magnetic field H. In analogy with the A^^^ case^'^, we expect 
various phases in the H — A plane, where A = ch 77. In particular, for A > 1 there should 
be a range of > for which the model exhibits massless behavior. 

We focus on the point A — > 1, if — > in the massless phase, where the model becomes 
SM(3)-invariant^^'^^. We solve the TBA equations in a systematic low-T and small-H 
perturbative expansion, along the lines of Johnson and McCoy^. These calculations employ 
Wiener-Hopf techniques, which were first used in a similar context by Yang and Yang^°. 
However, here we deal with a system of integral equations which requires factorization of a 
matrix kernel^^ . We then compute the free energy F(T, H), and determine the specific heat 
and magnetic susceptibility. Using the well-known relation^^'^^ between low-temperature 
specific heat and central charge (c) for critical systems, we arrive at the value c = 2 
for the SM(3)-invariant chain. This value coincides with that obtained^^ from finite-size 
corrections, and is expected from the equivalence of this model to the level-one su{3) WZW 
model in the continuum limit. We believe that our value for the magnetic susceptibility is 
new. 

While there are other ways^'^ of calculating the low-temperature specific heat within 
the general TBA approach, the method pursued here has the merit of treating this calcu- 
lation in the same manner as the one for the magnetic susceptibility. 
* in the sense defined above 
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We see no difficulty in computing thermodynamic quantities within the TBA approach 
at other points in the H — A plane. However, we are primarily interested in the massless 
phase, as this is where the connection with field theory is better understood. 

The outline of our paper is as follows. In Section 2, we formulate the TBA equations 

(2) 

for the closed A2 chain in the noncritical regime, with rj real. In Section 3, we take the 
limit 7/ — > 0, and arrive at the TBA equations for the antiferromagnetic stt(3)-invariant 
chain in the fundamental representation. In Section 4, we solve these equations for small 
values of T and H, and in Section 5 we calculate the free energy, specific heat, and magnetic 
susceptibility. We present further discussion of our results in Section 6. 



2. TBA equations for the A^^ chain 



(2) 

The Hamiltonian for the closed A^-site A2 chain is given implicitly in terms of the 
A'2^ R matrix in Eq. (1.1). (Explicit expressions are given in Refs. 15, 16, and 25.) In an 



external magnetic field H, the corresponding energy eigenvalues are 
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M 



E 



sh r\ 



H{N -M) 



(2.1) 



smr]{Xk - f)sinr7(Afe + |) 

where the real part of the complex variables have values in the interval [—iv/2r] , 7r/2?7] 
and satisfy the Bethe Ansatz (BA) equations 



sin?7(Afc + |) 
sin?7(Afe - f) 



TV 



M 



=-n 



smr]{Xk - Xj + i) cosr]{Xk - Xj - |) 
smr]{Xk - Xj - i) cosr]{Xk - Xj + |) 

k= 1,---,M. 



(2.2) 



We consider the (noncritical) regime where rj is real and (without loss of generality) posi- 
tive. 

We shall investigate the thermodynamics of this model, so we shall need to solve the 
BA equations in the A?" — > 00 and M — > 00 limit, with M/N fixed. As is customary, we 
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adopt the "string hypothesis" which states that all the solutions {A^ ,k = 1, • • • , 00} are 
collections of "strings" of "length" n of the form 



n + 1 



-I 



(2.3) 



where I = l,---,n; a — 0, 1,---,M„; n = l,---,oo; and the "centers" are real. A 
particular solution of the BA equations corresponds to a set of non- negative integers {M^} 
and the real numbers AJJ for each n. Observe that the total number of A variables, 
and BA equations, is M = Yl^=i nM^. These string solutions are^^ the same as those of 
the noncritical A^"^ chain^. 

It will prove convenient to introduce the functions 

sinr7(A+ 'f] 



Pn{X) = iln 
qn{X) = iln 



sin?7(A — 

cos ?7(A + 



vn \ 
2 ' 

in ^ 
2 ' 



cos?7(A - f ) 



(2.4) 



and the matrices 



min(n,m) — l min{n,m) 

=nm(A) =Pn+m(A)+2 ^ P\n-m\+2l{>^)+P\n-m\W+ ^ q2l-n-m-l{X) ■ (2.5) 

1=1 1=1 

We now substitute (2.3) into the BA equations (2.2), and take the product of the resulting 
equations for A^"^'^'' over the n values of Taking the logarithm, we then obtain the 
following equations for the centers X^: 



h-ix-) = j-, 



a = 0,l,---,M„; 



n = 1,2, 



(2.6) 



where 



( oo Mm 1 

^(A) = - ^ Npr,{X) -J2Y.^nm{X- Xf) \ 
[ m=l/3=0 J 



(2.7) 



and are integers or half-integers. We make the conventional assumption that h"'{X) is 
a monotonic increasing function of A. Let J denote the set of allowed values of J^, and 
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J its complement. If A is such that h'^{X) e J, it is said to correspond to a particle (of 
rapidity A). If A is such that h^{X) G J, it is said to correspond to a hole. Let Pn{X) be 
the density of particles and Pn(A) be the density of holes. Then 

Pn(A) +Pn(A) = lirn^ ^ • (2-8) 

An expression for the right hand side of this equation can be found by the substitution of 
iV"^ J2i3 by / d\'p{X') in (2.7). This leads to the equation 

oo 

+ On , (2.9) 



where 



and 



^nm('^) — 5nm^(A) + (1 ^nm) '^In— ml (A) 

min{n,m) — l 

(A) , (2.10) 

1=1 

min{n,m) 

Bnm{X)= b2l-n-m-l{X) , (2.11) 

1=1 



a„(A) = ^ Ap„(A) = 1 f'"") . (2.12) 
2iT dX TT ch(?7n) — cos(2?7A) 

, / X X 1 , , , ?7 sh(?7n) 

27r ctA TT ch(?7n) + cos(2?7A) 



and * denotes a convolution, 

»7r/2r7 
T/2r, 

For future reference we note here that 



/7r/2r7 
dX' f{X-X')g{X'). (2.14) 
-7r/2r7 



5nm(A) = -(s*Anm)(A+;^), (2.15) 

2r/ 



where s(A) is defined by 



s(A) = ^ y e-^'^'^—^. (2.16) 
^ ^ TT ^ 2ch(r)k) ^ ' 



We also note that s(A), which can be expressed in terms of the Jacobian elliptic function 
dn, has the property that 

s*a2 = ai-s, s * (an+i + an-i) = an n>l. (2-17) 

These and other relations which we give below can be easily derived with the help of 
Fourier transforms, for which we use the following conventions 

r>7I"/2?7 

e-'^^'Vfc , fk= I 

k=—oo 

The thermodynamic limit of the energy per site is 



/(A) = ^ V e-^'^^^h , h= dX e^'^^^fiX) . (2.18) 



while the entropy per site is 

/ dX [{pn + pn) In (pn + Pn) " Pn In Inp,] . (2.20) 



r-7r/2r? 



(2.19) 



AT 



The equilibrium value of p„ at temperature T is determined^ by extremizing the free energy 
per site F/N — {E — TS)/N. Note that the variation of pn is determined in terms of the 
set of variations {5pn} by the constraint (2.9) which implies that 

SpniX) ^ ^^^^^ _ ^ ^^^^^ _ ^2.21) 



Spm{y) 

Using this one finds that F/N is extremized when the functions 



satisfy the thermodynamic Bethe Ansatz (TBA) equations 

T In (l + e^"/^) = J2 i^nm + Bnm) * T In (l + e-^-"/^) - 27r^an + nH . (2.23) 

m=l ^ 

Using the TBA equations one finds that, in equihbrium, 

- = -T V / d\ an(A) In ( 1 + e-^"^^)/^ ) - i? . (2.24) 

As in other integrable models, the free energy can be re-expressed as a functional of only 
one of the en(A)'s. To see this, we introduce the matrix function 

KilrS}^ = <^(A)(^nm " s(A) (5n,m+l + (^n,m-l) , (2.25) 

where s(A) was defined in (2.16). As the notation suggests, it has the property 

oo 

{Kl' * ^n'm) (A) = ^{\)Km , (2.26) 

n' = l 

which follows from (2.17). It has the further properties 

oo oo 

^ {A-l * a J, (A) = s(A)5„i , 5] * m = . (2.27) 

m=l vn=\ 

Inserting A~^A into (2.24) and using (2.27) and (2.23) with Bnm replaced by the right 
hand side of (2.15), one finds that 

— = -27r^ / dX ai(A)r(A) -T dX r(A) In ( 1 + e^i^/^ ) . (2.28) 

N V J-n/2ri J-7r/2r, ^ ^ 

Here 

oo 

fc= — OO 

where Sfc = l/2ch.r]k are the Fourier coefficients of s(A). 

From (2.28) we see that in order to compute thermodynamic quantities of our model, 
we need only determine ei(A). Unfortunately, the TBA equations (2.23) are coupled non- 
linear equations for all the en(A). However, at low temperature these equations linearize 
and a determination of ei (A) becomes possible. 
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The type of low-temperature expansion depends crucially on the values of the anisotropy 
parameter rj and the magnetic field H, as shown for the A\^^ model by Johnson and 
McCoy^. We shall not pursue here a similar exhaustive analysis of the model. In- 
stead, we shall concentrate on the limits 

(2) 

In the limit 77 — > 0, the A2 model reduces to the sw(3)-invariant quantum spin chain in 
the fundamental representation^^'^^. In the following section, we shall explain how the 

(2^ 

TBA equations of the SM(3)-invariant model emerge from those of the A2 model. We 
remark that the su(3)-invariant model allows the introduction of two external (magnetic) 
fields because su{3) has rank two. By our limiting process we find only one combination. 



3. The r]^0 limit 



We recall the A^^^ BA equations (2.2): 



sinr7(Afc + f) 
sinr/(Afe - |) 



N 



M 



=-n 



sin?7(A;s - Xj + i) cosr]{\k - Xj - |) 
^^j^ smr]{Xk - Xj - i) cosr]{Xk - Xj + |) 

k= 1,---,M. 



(3.1) 



Consider a solution {Ai,- • •, Am}- For 77 — > 0, some of the A^'s remain finite. We call 
these solutions X^^\ k = 1, ■ ■ ■ , M^^\ The remaining solutions become infinite; of these, 
we restrict ourselves (following Ref. 16) only to those solutions whose behavior as 77 — > 
is given by ±^ -|- X^^\ k = 1, • • • , M^'^\ with A^^-* finite. Rewriting (3.1) in terms of the 
new variables A^^"* and A^'^'*, we obtain two families of BA equations, which become in the 
?7 — > limit 



,(1) 



N 



X 



(1) 



Wl) , M(2) x(l) 3,(2) , 

^(l) _ ^(1) .. 11 ^(l 



.1 A 
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(1) 



x(2) , i 
Xj, +2 



k 



M^^^ \(2) _ \W _ i M(2) ,(2) _ x(2) . 

11.(2) ,(1) , 11 a(2) x(2) .' '^^^ ^^-^^ 

These are precisely the BA equations for the stt(3)-invariant chain^^. 

In the formulas of the previous section, we must now distinguish two classes of string 
centers. The string hypothesis becomes 



where / = 1, • ■ • , n; a = 0, 1, ■ ■ ■ , M^^; r = 1, 2; n = 1, ■ ■ ■ , oo; and the centers A^'^'^'* are 
real. The total number of A's of type r is given by 

oo 

MW = ^nMW, r = l,2. (3.4) 

n=l 

We correspondingly have densities of particles pn\X) and holes Pn\x)- The integral 
equations satisfied by these densities are 

oo 2 

Pn^ * * = anSrl , (3.5) 

m=l s=l 

where Anm and a„ are the corresponding ry — > limits of (2.10) and (2.12), 



1 poo / I I 



oo 

oo 



^„m(A) = ^/ <jwe-""lcthJ|i| e-^l"-'"l-e-'*("+'") , (3.6) 



= ilL*^ 

C^s are the components of the 2x2 matrix 



where 



1 /""^ 1 1 

s(A) = ^ dw e"*^'^— — = — -— , (3.9) 
^ ' 27r 2chf 2ch7rA ' ^ ' 
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and * now denotes the convolution 



/oo 
dX' f{X-X')g{X'). (3.10) 
-oo 



The TBA equations are 

oo 2 

rin(^l + e^-'/^) = J2 ^^nm*a. *7^1n(l + e-^- - 27ra„5^i + ni? , (3.11) 



m=l s=l 



where 

e(:HX) = Tln(^] . (3.12) 

\Pn (A)/ 

Forming the convolution of the TBA equations with (which is given by the expression 
(2.25), with s(A) now given by (3.9)) and with C~J', we obtain 

oo 2 

Tln(l + e-^'VT) ^Y.Y.^nL*C-UTln{l + e^'^'/^)+2n6r,,s(^\ (3.13) 

m=l s=l 

where 

^ ^ ^ ^ 27ry_<^ 2chf ^1 V3 sh27rA ^ ^ 

The equilibrium free energy is given by 

771 ^ pOO 

I^ = -^E/ rfAa„(A)ln(l + e-^^'(^)/^)-iy. (3.15) 
With the help of (3.13), this expression can be cast in a form which depends only on 



F 
N 



/oo 2 „oo 

dXai{X)s^^\X)-Tj2 rfA s('-)(A)ln(l + e^i ^^)/^j . (3.16) 
-oo ^— oo 



In the following sections, we shall solve the TBA equations (3.11) for small values of T 
and H, and evaluate the expression (3.16) for the free energy. 
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4. T and H expansion 

We begin by rewriting the TBA equations of the su (3) -invariant chain (3.11) forn = 1 
in the form 

2 oo 2 

Tln^l + e^i"^'/^) =^Aii*a,*rin (1 + 6"^^^'/^) + J] J]Ai^*a,*Tln (^l + e"^-/^) 

s=l m=2 s=l 

- 27rai5^i + ii" . (4.1) 

We make the crucial assumption that e"m {\) > for m > 1. (The results which we obtain 
below do not contradict this assumption. For other examples, see Ref. 9.) This means 
that exp(— em /T) goes to zero exponentially as T — > for m > 1 and can therefore be 
neglected. On the other hand, e^'^''(A) can have either sign. Defining £('^)(A) to be the 
T^O limit of e^{\\), 



we see that 



£W(A) = ^im er(A), (4.2) 



hrn^ Tln(l + e^'i"'/^) = ±6^'^)^ , (4.3) 
where we use the standard notation 

£~ = ^(£-|£|), £+=£-£". (4.4) 

It follows from (4.1) that 

2 

g(r)+ = _ ^ ^ Crs * e^'^- - 2TTai5rl + H . (4.5) 
s=l 

It will be convenient for subsequent analysis to work with the following equivalent equa- 
tions, which do not involve £*^'')~: 

£(1) = H- 27rs(i) + h * £^1)+ + g * e^^^+ , 

£(2) = H- 27rs(2) ^ ^ ^ ^(2)+ ^ g ^ ^(1)+ ^ (4_g) 
13 



where s^^^ and s^^^ are defined in (3.14), and 
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-sW + si'^) ^ ai , /i = + s« * ai . (4.7) 



In order to obtain (4.6), we substitute into (4.5) the exphcit expressions for An and Crs 
as given by Eqs. (2.10) and (3.8), as well as = e*^^) — and then we solve for s^'^'' 

in terms of £('^^+ with the help of Fourier transforms. 

We observe for future reference that by a similar procedure, (4.1) can be recast (after 
neglecting terms involving with m > 1, but before taking the T — > limit of ei^^) as 
follows: 

e(') = H- 27rs(i) + h*T\n{l + e^^'V^^ + ^ * Tin (l + e^i'^^^ ^ 

ei') = H- 27rs(2) + /i * Tin (l + e^i'^T^ + ^ * Tin (l + e^i'^T^ _ (43) 

Evidently, the T — > limit of these equations gives (4.6). 

We first briefly consider the case H — 0, which corresponds to the ground (vacuum) 
state. Since s^^\X) and s^'^\X) are positive for all A, the equations (4.6) have the solution 

eW(A) = -2..W(A) = -^^^i^i^, 
^ ^/3 ch(7rA) 

e(^)(A) = -27r.(^)(A) = • (4.9) 

V3 sh(7rA) 

From the definition (3.12) and the constraint equations (3.5), and also from the assumption 

(r) 

that Sm > for m > 1, it follows that the ground-state densities of particles and holes are 
given by 

(1) 1 ch(7rA/3) 1 M/3) 

pW(A) = 0. (4.10) 
14 



That is, the ground state of the antiferromagnetic S'u(3)-invariant chain consists of a "con- 
densate" of strings of length 1. Defining momenta p^'^\X) of the quasi-particles by (see, 
e.g., Kirillov and Reshetikhin^^) 

^p(^)(A) = eW(A), (4.11) 
we see from (4.9) that for A — > oo, there is a hnear dispersion relation 

= vsp^'''^ , (4.12) 

with the velocity of sound Vg = 27r/3. 

We now turn to the case H 0. For H 0, Eqs. (4.6) can in principle be solved by 
iteration, with the zeroth-order solution given by 

e''o\X) = H-2TTs^''\X). (4.13) 

Since the functions s^'^\X) are positive and monotonically decreasing, for — > the 
functions e^'^\X) have a single zero. This zero is for A = O (InH). Assuming an expansion 
in powers of {liiH)~^, we conclude that 

£W(aW) = (4.14) 

for 



(r) 3 
271 



2.-y -vini^,' ' ^'-''^ 

where the constants k^'^'' (which are independent of H) have still to be determined. We 
shall further assume that the functions ^'■'^''(A) have no other zeros in the interval (0 , oo). 
Hence, £*^'')(A) < for A in the interval (0,q;*^^^), and £('^)(A) > for A in the interval 
(aW ,oo). 

Observe that Eqs. (4.6) are temperature-independent. Using the same approximations 
in the expression (3.16) for the free energy this too will be temperature-independent. To 
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find the leading order temperature dependence we need to compute the leading correction 
to the solutions e^*^"* = e^*"^ of the linearized equations (4.6). In order to obtain this 
correction, we make the substitution 

eM(A) =£M(A)+?7^'^)(A) (4.16) 

in (4.8) and expand to leading order in rj^'^K Since e^"^^ is a solution of the linearized 
equations (4.6), we shall find inhomogeneous terms in the resulting equations for r]^'^\ 
Indeed, we have that 

rj^'^=h*{Tln [l + e(^^''+^^'^)/^] - s^'^+] + g * {t In [l + e(^''^+^^'')/^] - e^^^+] , 

(4.17) 

and the similar equation for r]^"^^ is obtained by interchanging the superscripts (1) and (2). 
Because £*^''^(A) is an even function of A with a single zero, at A = a^'^\ for positive A, we 
see that (assuming ij^'^^ /T is small and keeping terms linear in ij^'^^) 

/*{Tln [l + e(^'^'+^'^')/^] -eW+} 



-a 



(r) 



CXD 



+ / ^ JdA' /(A-A'){Tln[l + e(=^^^(^')+^''^'(^'))/^] -eW(A')} 



+ / dX' /(A- A')Tln 



(r) 

— ^ poo 
+ / 



1 + e 



(e('-)(A')+77(^>(A'))/T 



(4.18) 



dX' f{X-X')r]^'-\X')+Ef^ 

(r) 

where the inhomogeneous term is given by 

4'^)=/*Tln(l + e-l^^'^^l/^) . (4.19) 

(Here / represents either of the kernels h and g appearing in (4.17).) For T — > 0, the major 
contribution to the integral comes from the regions near the zeros of e^'^\ so we expand 
£W(A) about A = aW, 



£W(A) = (a - + O ((A - a(^))2) , 
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t(r) 



dX 



(4.20) 



One then finds that the leading T- dependence of e'J^ is 



E^p{X) = 



/(A-«(^)) + /(A + a('^))J j du ln(l + e-'^) 
/(A-aW) + /(A + aW)] . 



(4.21) 



Taking into account these 77 - independent terms in the expansion of (4.8) to 0(77), we 
obtain the following linear integral equation for ri^^\ 



r7W(A) 



+ 

+ 
+ 



+ / \d\' h{X-X')r^^^\X') 

-00 JaW J 

dX' g{X-X')r]^^\X') 



— ^ poo 

+ / 



6t(i) 
6t(2) 



/i(A-aW) + /,(A + «(!)) 
'g{X-a^^^)+g{X + a^^^) 



(4.22) 



The similar expression for ry^^^ is obtained by interchanging the superscripts (1) and (2) 
on ?7, a and t. These equations for 77^^) and those of (4.6) for e^*") complete our results for 

(r) 

the T - expansion of e^^ . We now turn to the expansion in powers of In H, by which from 
Eqs. (4.6) and (4.22) we shall generate systems of integral equations of the Wiener-Hopf 
type. 

It will prove convenient to work with the functions 



- t A < ' 



and 



r('^)(A) 



^-^^;^r,(r)^X + a(r)) A>0 



(4.23) 



(4.24) 



2 in 



A<0 

instead of the functions £(^)(A) and //^(A). The factors e^'^^^^'^s and (e'^^^^"^^ /^'^ /T' 
(4.23) and (4.24), respectively, are chosen such that the driving terms in the equations for 
5(^)(A) and T(^)(A) have a nonvanishing limit as T ^ and ^ 0. The factors of k^^), 
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which were first introduced in (4.15), appear in (4.23) for a reason which will be explained 
below. 

We return now to Eq. (4.6). We write the limits of integration explicitly, keeping 
in mind that ^('^^(A) > for — oo < A < —a^'^^ and for a^'") < A < oo; and we shift 
the integration variables so that they run from to oo. Observe now that h{X + 2a('^)) 
vanishes as H ^ («('") — > oo) for finite A, and similarly for g (A + a(i) + a(2)). The 
functions g [X + a*^^-* — ct'-^'') and g [X + a*-^'' — ct*-^-*) remain finite, however, and survive in 
the H ^ limit of the equations for S^'^\X), which are (for A > 0) 

5W(A) = ^fl-«We-2-AM 



a/3 V 



+ 



dX' 



J 



h{X - X')S^^\X') + g{X -X' + a^^) - a^^^)S^^\X') 



5(2)(A) = ^(^l_«(2)e-2^A/3^ 



V3 



+ 



dX' 



'g{X -X' + _ a(i))^(i) (A') + h{X - X')S^'^^ (A')] • (4.25) 



Similarly, from (4.22), we see that the — > limit of the equations for T*^^^(A) are 
(for A > 0) 



t(i)(a) 



h{X) ^(A + a(i) -a(2)) 



+ 



5(2)' (0) 



dX' 



h{X - A')T(i)(A') + ^(A - A' + - a(2))T(2)(A') 



T(2)(a) 



h{X) ^(A + a(2) _c,(i)) 



5(2)' (0) 



+ 



+ 



dX' 



g{X -X' + a(2) - a(i))T(i)(A') + h{X - A')T(2)(A')1 , (4.26) 



where 



5(0' (0) = -^5(^)(A) 
dX 



g2,raW/3^(r)^(r) 



A=0+ 



(4.27) 



These equations can be written in the standard Wiener-Hopf form 

5('') (A) = (A) + ( A) + ^ / dX' K"-' (A - A') S^'^ (A') 
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2 nCXD 

s=l"'-oo 

- oo < A < oo , (4.28) 



and similarly 



where the kernels K^^{X) are the components of the 2x2 matrix 

K(X)-( ^(A + aW-«(2))\ 

and 

^.-.w.jaKi-"';--^'^) ^>«, (,30) 

\ A<0 

/jf)(A) = <^ S(2)'(0) ^ S(i)'(0) A>U .^33. 

\ A<0 

,(r)/^s. _ / A > r/i 

W - I _ ^2^^ j-^ rfA' K-(A - A') T(^)(A') A < • ^^-"^^^ 

We shall solve these equations by Fourier transform. We define the Fourier coefficients 

of 5('")(A) and tW(A) by 

/oo l-OO 
dX e'^"" S^''\X) , f^''\u)= / dA e*^'^TW(A). (4.35) 

-00 J —00 

Since ;S('")(A) and T('")(A) vanish for A < 0, the functions ^('"^(a;) and f ('")(a;) are analytic 

in the upper-half-plane Im u >0, which we denote by n_|_. Observe that since S^'^\0) — 

(as follows from e^'^\a^'^^) = 0), we have by contour integration 

5^(0) = -i lim ujS^''\uj) = , (4.36) 

loll— > 00 

where the limit is taken in 11+ . We also note that 

S^^y (0) = (A) =- hm a;2^W {u) . (4.37) 

CIA A=0-|- a» — >oo 
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Care must be exercised in the derivation of this result because of the discontinuity of the 
derivative of S^'^'^ (A) at A = 0. 

The Wiener-Hopf equations for S^'^^ in Fourier space are 

2 



^^^n^) = /i^^M + ^i'^M + E^"H -^^'^M' (4.38) 

8 = 1 



where ^''^(a;) are the components of the 2x2 matrix 

Observe that (for uj and o:*^^) real) this matrix is Hermitian, 

K(a;)t = k{u) . (4.40) 

The factors of k*^^^ in the definition (4.23) of S^''^ were chosen to arrange for this to be the 
case. 

The Wiener-Hopf equations for T^"^^ in Fourier space are similarly found to be 

2 

TM(a;) = f!^\u^) + hP{u^) + ^^'^('^) • (4-41) 

Since ^1 — K{ijj)^ is nonsingular, Hermitian and positive-definite at a; = 0, it is 
positive-definite for —oo < u) < oo. Thus, Theorem 8.2 of Gohberg and Krein^^ implies 
that the following factorization exists 

{l - K{lo)^ ^ = G+icj) G-{u), -oo < a; < oo , (4.42) 

where G^{uj) and G^^{u}) are analytic in H+ with G^{u>) ^ 1 as a; — > oo in H+, and 
G-{ijj) and Gz}{u:) are analytic in H_ with G-{ijj) — > 1 as o; — > oo in H_. Moreover, the 
fact K{—\) = K{X)'^ implies that (for uj in H_) 

G'_(a;)^ = G'+(-a;), (4.43) 
20 



where the superscript T denotes transpose. This lemma can be proved from the formulas 
developed in Ref. 31. As we shall see, explicit expressions for and G- are not needed 
to compute the free energy to the order in which we work. A similar phenomenon occurs 
in the Wiener-Hopf calculations of Yang and Yang^° and Johnson and McCoy^. 

Using the factorization (4.42), the Wiener-Hopf equation (4.38) for 5'*^^^ can be rewrit- 
ten (in 2 X 2 matrix notation) as 

G^' S = G_(^fs + hs) . (4.44) 

Observe that the left hand side is analytic and bounded in n_|_, whereas G-bs is analytic 
and bounded in n_. The term G-fs has a decomposition as 

gJs = P- (G-fs) + P+ (gJs) , (4.45) 

where P± (^G-fs) is analytic in II-i-. This decomposition is uniquely specified by the 
requirement P± (G-Js) — > for a; — > oo in Il±. Taking the P+ projection of (4.44), we 
have that 

S = G+P+ (gJs) . (4.46) 

From (4.30) we compute that 

fPi.) = . (4.47) 

where one is to take e — > at the end. The decomposition (4.45) of G-fs is then found 
by subtracting the residues of fs, i.e., 

G.(.)M.) =^ {-1- - G.(-ie)) (;) - (G.M - G-(-2.V3)) (^S:! )} 

- Ts l^'^-'-' ( i ) - ;:t^«-(--/^' ) } ■ 

Hence, 

(4.49) 
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The boundary condition (4.36) can now be seen to be equivalent to the condition 



(-271^/3) (^S) =G_(0)Q) , (4.50) 

which can be solved for the parameters k^'^\ Using this result we conclude that Siuj) is 
given by 

S{uj)^'^(—^. ]-—\g+{u)G-{Q)(\\. (4.51) 

From Eq. (4.37) and the fact that G+(a;) — > 1 as |a;| — > oo in 11+ , we find that 

^^'"^'(0) = ^G'_(0)Q) • ^^-^^^ 

We turn now to the equation (4.41) for T{u)). Proceeding as before, we use the 
factorization (4.42) to arrive at the formal solution 

f = G+P+ (g_/t) (4.53) 

(cf. Eq. (4.46) ). The explicit calculation of f^'' is difficult, but can be avoided by the 
following trick. Consider the functions f^^\X) defined as 

fW(X) = ^(^) + ^(A + aW-«(^)) 
^ SW{0) 5(2)' (0) 

^ ^ ^ >S(2)'(0) ^ 5(1)' (0) ' ^ ^ 

where A ranges over the entire real line. From Eqs. (4.32) and (4.33), it is evident that 

f!f\x) = ft\x), (4.55) 

where 

(4.56) 
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The Fourier transform of f^'^\X) is readily computed, and can be expressed in terms of 
the kernel K, 

m = k{u) (l^^S'foj) . (4.57) 
From the factorization (4.42), it follows that 

/++/_=(i-G:'G-)(;/fs: {;;>). (4.58) 

After multiplying both sides of this equation by G_, we see that the P+ projection of 

G-f+ is given by 

p+(G_A)=a-G-(;/^|j;:[°]) , (4.59) 

where a is a constant. Requiring the right hand side to vanish for |a;| oo in 11+ 
determines this constant to be 

We conclude from (4.53) that T{u>) is given by 

r(a;) = (G+(a;)-l)(^J/;^j2r[o|) , (4-61) 

where S'^^)'(O) is given by (4.52). 

To summarize this section: we have made the expansion (4.16) of e^'^'*(A), and we have 
changed in Eqs. (4.23), (4.24) from the variables e^''\X), rj^'^^X) to the variables <S(^)(A), 
T('")(A), respectively. Using Wiener-Hopf methods, we have determined in Eqs. (4.51), 
(4.61) the corresponding Fourier transforms S^'^\uj), T^'^\u>) in the limits T — > and 
— > 0. These expressions involve G+{(jj) and G_(a;), which appear in the factorization 
(4.42). In the next section, we shall use these results to calculate the free energy. 
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5. The free energy 



Substituting the expansion (4.16) of e^*^-* into the expression (3.16) for the free energy 
(keeping in mind the discussion immediately following (4.16)), we obtain 



2^2 



N 



= eo 



E 



3 ^ t(^) 

r=l 



r=l 



2J2 rfAs(^)(A) p^)(A) + ?7(^)(A)] , (5.1) 



where the ground state energy per site eo is given by 

/oo 
-oo 



(5.2) 



We aim for a double expansion of F/N to quadratic order in both T and H. To this 
end, we rewrite (5.1) in terms of S^'^\X) and T^'^\X). Consider first the term. For 
— > 0, we can make the approximation 



1 



and hence this term becomes 



E 



3V3 fr[ s(^y (0) ■ 

We next manipulate the integral in the last term as follows 

/■oo 

/ dXs^''\X) £('^)(A)+?7('^)(A) 
dX s^'^^A + a^")) 



(5.3) 



(5.4) 



-27ra('^>/3 



(r) 



„27^2^(r) 



(6.5) 



gg-27ra('-)/3' 

Taking the a^'^^ — > oo limit of this expression, we arrive at the following expression for 
F/N, 

x/3if2 2 



F 
N 



T7 = eo 



27r2 



r=l 



^^g-2.A/3^(r)(;^) 



(0) 



(5.6) 
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Writing the functions S^'^\X) and T(^)(A) as the Fourier transforms of ^(a;) and ^(0;), 
respectively, and then performing both the A and u> integrals, we obtain 



Co A —B 



N 



27r2 



(5.7) 



with 



r=l 
2 



K 



(r) 



= K 



T(27rz/3)+K^^(2)'^0) 



(5.8) 



where we again switch to a matrix notation. 

Since both A and B involve , we begin by observing from (4.50) that 



''=(^S)=G_(-27rV3)-^G_(0)(^}) . 



(5.9) 



Taking the transpose of this equation, and using the property (4.43), we obtain 



K, 



= («W «(2))=(1 1) G+(0)G+(27rz/3)-^ 



(5.10) 



We now consider A. Evaluating S{27ri/3) using the expression (4.51), we obtain 

A=y^{l 1) G+(0)G_(0)Q) . (5.11) 

Recalling the factorization equation (4.42) and the explicit expression (4.39) for K{u;), we 
conclude that 



V3 



(1 1) (i-K(o)j (l)=Vs 



(5.12) 



There remains to compute B. Evaluating T(27rz/3) using the expression (4.61), we 
obtain (after a crucial cancelation) 



B=iii) G.(o)(;/^!:y[°]) 



(5.13) 
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Recall Eq. (4.52), 

Taking the transpose of this equation and again using the property (4.43), we see that 



B 



3\/3 
4^ 



.-(0) .-(0)) (y^;^^^ 



(5.15) 



The expression for the free energy per site is therefore 



' rj-i2, 



(5.16) 



As foreseen above, this result was obtained without using explicit expressions for G+ and 
G-. It follows that the magnetic susceptibility and specific heat, to lowest order, are given 

by 



X = 



Ch = -T 



df^ [n 



T TT^ 



H 



(5.17) 



For a critical chain, the free energy per site is given by^^'^^ 



TTC 



(5.18) 



where c is the central charge and Vg is the velocity of sound. For our model, Vs — 2tt/3 
(see Eq. (4.12)), and hence c = 2. 
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6. Discussion 



For an integrable model, the algebraic calculations leading to the energy eigenvalues 
in terms of solutions of the BA equations are quite elegant and precise. The same can- 
not be said for the corresponding thermodynamic calculations, at least in their present 
formulation. This is both surprising and disappointing. For instance, one expects that 
there should be an analogue of the Sugawara construction for integrable models, by which 
one could determine the low-temperature specific heat (central charge) by purely algebraic 
means. The fact that an explicit Wiener- Hopf factorization is not needed to compute such 
properties also suggests that an alternative approach should be possible. Indeed, very re- 
cently, progress has been made^^ towards an algebraic formulation of the thermodynamics 
of integrable models. 

Despite its short-comings, the approach which we have followed here to investigate the 
low-temperature thermodynamics of the closed stt(3)-invariant chain in the fundamental 
representation is nevertheless practical. It reproduces the known result for the central 
charge, and evidently, it can be implemented for su{n). 

We recall^^ that, in the continuum limit, integrable spin-s sM(2)-invariant chains"^^ 
are described by level k = 2s su{2) WZW models. Moreover, Affleck has found^^ a 
simple relation between the magnetic susceptibility x ^md the level, namely x'^s = k/2TT. 
Presumably, there is a generalization of this relation to the su{n) case, for which there are 
n — 1(= rank of su{n)) magnetic susceptibilities x^, z = l,---,n — 1. We have determined 
for sm(3) a particular linear combination of magnetic susceptibilities, which is dictated by 
the special imbedding so(3) C su{2>) which characterizes the A)^ chain. Our result may 
be relevant to the stt(n)-generalization of Affleck's relation. 

(2) 

Having treated the noncritical regime of the A2 chain, the task now is to investigate 
the critical regime. As noted in the Introduction, this may be more feasible for the open 
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chain, which has Uq[su{2)] symmetry. We hope to report on this problem in the future. 

We thank A. Jacob for bringing Ref. 31 to our attention, and E. Melzer for useful 
discussions. Part of this work was performed at the Aspen Center for Physics. This work 
was supported in part by the National Science Foundation under Grants No. PHY-90 
07517 and PHY-92 09978. 
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